
log using "Log_replication_results.log", replace

********************************************************************************
* Section "EXPERIMENTAL DESIGN", subsection "Procedures"						
********************************************************************************
	
*	Information on participants	
		quietly do "Data_Reading_subjects.do"

			sum age if period==2
			* 24.2
			
			tab sex if period==2
			* 94 female / 86 male
			
			tab student if period==2
			* 177 students
			
			tab econ_student if period==2
			* 77/180 (42.78%) econ students
			
			sum gesamtauszahlung if period==2
			* 15.17 EURO
		

********************************************************************************
** Section "EXPERIMENTAL RESULTS" 						                
********************************************************************************
	
*	Length of committee decisions
		quietly do "Data_Reading_decision_dynamics.do"		
				
		* Duration of group decisions		
				sum time_in_stage2 if decision==4
				list time_in_stage2 if decision==4
		
		* Number of proposers in group decisions
				sum num_ag if decision==4			
	
	
* 	Did committees choose points inside the uncovered set?
		
		* UCS has fairness properties (Setting A)
			sum final_ucset if decision==4 & (treatment==1 | period==1)
		* UCS lacks fairness properties (Setting B & C)
			sum final_ucset if decision==4 & (treatment!=1 & period==2)
			

		// TREATMENT 1, Period 1
			tab final_ucset if treatment==1 & period==1 & decision==4

			/* 
			Decisions inside UC-set: 11 , Decisions outside UC-set: 1 
			Size uncovered set: 313 units
			Size Pareto set: 780 units
			percentage of outcomes occurring in the uncovered set: 91.67%
			percentage of outcomes that one would expect to find in the uncovered set 
			given a uniform distribution of outcomes throughout the Pareto set: 40.13%
			*/	

			bitesti 12 11 .40128205
			
		
		// TREATMENT 1, Period 2
			tab final_ucset if treatment==1 & period==2 & decision==4

			/* 
			Decisions inside UC-set: 12 , Decisions outside UC-set: 0 
			Size uncovered set: 313 units
			Size Pareto set: 780 units
			percentage of outcomes occurring in the uncovered set: 100%
			percentage of outcomes that one would expect to find in the uncovered set 
			given a uniform distribution of outcomes throughout the Pareto set: 40.13%
			*/	

			bitesti 12 12 .40128205
			
			
		// TREATMENT 2, Period 1
			tab final_ucset if treatment==2 & period==1 & decision==4

			/* 
			Decisions inside UC-set: 11 , Decisions outside UC-set: 1 
			Size uncovered set: 313 units
			Size Pareto set: 780 units
			percentage of outcomes occurring in the uncovered set: 91.67%
			percentage of outcomes that one would expect to find in the uncovered set 
			given a uniform distribution of outcomes throughout the Pareto set: 40.13%
			*/	

			bitesti 12 11 .40128205
			
			
		// TREATMENT 2, Period 2
			tab final_ucset if treatment==2 & period==2 & decision==4

			/* 
			Decisions inside UC-set: 1 , Decisions outside UC-set: 11 
			Size uncovered set: 313 units
			Size Pareto set: 780 units
			percentage of outcomes occurring in the uncovered set: 8.33%
			percentage of outcomes that one would expect to find in the uncovered set 
			given a uniform distribution of outcomes throughout the Pareto set: 40.13%
			*/	

			bitesti 12 1 .40128205
			
			
		// TREATMENT 3, Period 1
			tab final_ucset if treatment==3 & period==1 & decision==4

			/* 
			Decisions inside UC-set: 11 , Decisions outside UC-set: 1 
			Size uncovered set: 313 units
			Size Pareto set: 780 units
			percentage of outcomes occurring in the uncovered set: 91.67%
			percentage of outcomes that one would expect to find in the uncovered set 
			given a uniform distribution of outcomes throughout the Pareto set: 40.13%
			*/	

			bitesti 12 11 .40128205
					
			
		// TREATMENT 3, Period 2
			tab final_ucset if treatment==3 & period==2 & decision==4

			/* 
			Decisions inside UC-set: 0 , Decisions outside UC-set: 12 
			Size uncovered set: 313 units
			Size Pareto set: 780 units
			percentage of outcomes occurring in the uncovered set: 0%
			percentage of outcomes that one would expect to find in the uncovered set 
			given a uniform distribution of outcomes throughout the Pareto set: 40.13%
			*/	

			bitesti 12 0 .40128205
				
	
* Figure 4. Spatial distribution of final results

			*Treatment 1, Period 1
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter final_y final_x if treatment==1 & period==1 & decision==4 & final_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter final_y final_x if treatment==1 & period==1 & decision==4 & final_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==1, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Decisions inside Uncovered Set"  35 "Decisions outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 1, Period 1}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("11 decisions inside Uncovered Set" "Average decision: (48.3|40.4)" "Standard distance deviation: 3.55", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment1P1, replace)

			*Treatment 1, Period 2
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter final_y final_x if treatment==1 & period==2 & decision==4 & final_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter final_y final_x if treatment==1 & period==2 & decision==4 & final_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==1, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Decisions inside Uncovered Set"  35 "Decisions outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 1, Period 2}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("12 decisions inside Uncovered Set" "Average decision: (49.7|41.0)" "Standard distance deviation: 2.29", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment1P2, replace)

			*Treatment 2, Period 1
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter final_y final_x if treatment==2 & period==1 & decision==4 & final_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter final_y final_x if treatment==2 & period==1 & decision==4 & final_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==2, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Decisions inside Uncovered Set"  35 "Decisions outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 2, Period 1}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("11 decisions inside Uncovered Set" "Average decision: (49.4|40.6)" "Standard distance deviation: 5.89", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment2P1, replace)

			*Treatment 2, Period 2
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter final_y final_x if treatment==2 & period==2 & decision==4 & final_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter final_y final_x if treatment==2 & period==2 & decision==4 & final_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==2, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Decisions inside Uncovered Set"  35 "Decisions outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 2, Period 2}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("1 decision inside Uncovered Set" "Average decision: (62.6|41.2)" "Standard distance deviation: 10.86", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment2P2, replace)

			*Treatment 3, Period 1
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter final_y final_x if treatment==3 & period==1 & decision==4 & final_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter final_y final_x if treatment==3 & period==1 & decision==4 & final_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==3, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Decisions inside Uncovered Set"  35 "Decisions outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 3, Period 1}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(small)) ///
						note("11 decisions inside Uncovered Set" "Average decision: (48.8|39.5)" "Standard distance deviation: 4.23", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment3P1, replace)

			*Treatment 3, Period 2
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter final_y final_x if treatment==3 & period==2 & decision==4 & final_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter final_y final_x if treatment==3 & period==2 & decision==4 & final_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==3, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Decisions inside Uncovered Set"  35 "Decisions outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 3, Period 2}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("0 decisions inside Uncovered Set" "Average decision: (32.0|41.9)" "Standard distance deviation: 3.43", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment3P2, replace)

			
			* COMBINE ALL GRAPHS	
				grc1leg Treatment1P1 Treatment1P2 Treatment2P1 Treatment2P2 Treatment3P1 Treatment3P2, ///
				graphregion(color(white) margin(vsmall)) ///
				graphregion(color(white) margin(vsmall)) ///
				cols(2) ///
				name(Figure_4, replace)
				graph display, ysize(14) xsize(7.8) scale(1) 

				graph close Treatment1P1 Treatment1P2 Treatment2P1 Treatment2P2 Treatment3P1 Treatment3P2 

				
				
* Spatial distribution of proposals				

		bys treatment period: sum prop_x prop_y if proposer_dec==1
		
		* Total number of proposals by treatment
		bys treatment period: count if proposer_dec==1
		
		* Number of proposals in UCS by treatment
		bys treatment period: count if proposer_dec==1 & prop_ucset==1
		
		* Number of proposals to the right of UCS in Treatment 2
		count if proposer_dec==1 & prop_ucset==0 & treatment==2 & period==2
		count if proposer_dec==1 & prop_ucset==0 & prop_x>50 & treatment==2 & period==2
		
		* Number of proposals to the left of UCS in Treatment 3
		count if proposer_dec==1 & prop_ucset==0 & treatment==3 & period==2
		count if proposer_dec==1 & prop_ucset==0 & prop_x<50 & treatment==3	& period==2	
		
		
		

* Figure 5. Spatial distribution of proposals

			*Treatment 1, Period 1
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter prop_y prop_x if treatment==1 & period==1 & proposer_dec==1 & prop_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter prop_y prop_x if treatment==1 & period==1 & proposer_dec==1 & prop_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==1, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Proposal inside Uncovered Set"  35 "Proposal outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 1, Period 1}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("17 of 28 (60.7%) proposals inside" "Uncovered Set" "Average proposal: (50.9|44.8)", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment1P1, replace)

			*Treatment 1, Period 2
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter prop_y prop_x if treatment==1 & period==2 & proposer_dec==1 & prop_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter prop_y prop_x if treatment==1 & period==2 & proposer_dec==1 & prop_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==1, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Proposal inside Uncovered Set"  35 "Proposal outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 1, Period 2}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("25 of 29 (86.2%) proposals inside" "Uncovered Set" "Average proposal: (49.4|41.8)", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment1P2, replace)

			*Treatment 2, Period 1
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter prop_y prop_x if treatment==2 & period==1 & proposer_dec==1 & prop_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter prop_y prop_x if treatment==2 & period==1 & proposer_dec==1 & prop_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==2, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Proposal inside Uncovered Set"  35 "Proposal outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 2, Period 1}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("26 of 36 (72.2%) proposals inside" "Uncovered Set" "Average proposal: (48.9|40.6)", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment2P1, replace)

			*Treatment 2, Period 2
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter prop_y prop_x if treatment==2 & period==2 & proposer_dec==1 & prop_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter prop_y prop_x if treatment==2 & period==2 & proposer_dec==1 & prop_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==2, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Proposal inside Uncovered Set"  35 "Proposal outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 2, Period 2}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("3 of 35 (8.6%) proposals inside" "Uncovered Set" "Average proposal: (63.3|39.4)", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment2P2, replace)				
				
			*Treatment 3, Period 1
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter prop_y prop_x if treatment==3 & period==1 & proposer_dec==1 & prop_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter prop_y prop_x if treatment==3 & period==1 & proposer_dec==1 & prop_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==3, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Proposal inside Uncovered Set"  35 "Proposal outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 3, Period 1}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(small)) ///
						note("20 of 39 (51.3%) proposals inside" "Uncovered Set" "Average proposal: (55.3|46.2)", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment3P1, replace)

			*Treatment 3, Period 2
				twoway ///
					(function y=51, range(49 49) color(gs12) clwidth(.8) aspect(1)) ///
					(function y=50.5, range(46.5 51.5) color(gs12) clwidth(.8)) ///
					(function y=50, range(44 54) color(gs12) clwidth(.8)) ///
					(function y=49.5, range(41.5 56.5) color(gs12) clwidth(.8)) ///
					(function y=49, range(39 59.5) color(gs12) clwidth(.8)) ///
					(function y=48.5, range(38 60) color(gs12) clwidth(.8)) ///
					(function y=48, range(37.5 60.3) color(gs12) clwidth(.8)) ///
					(function y=47.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=47, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=46, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=45, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=44, range(37.5 60.2) color(gs12) clwidth(.8)) ///
					(function y=43.5, range(37.5 60.8) color(gs12) clwidth(.8)) ///
					(function y=43, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=42, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=41, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40.5, range(38 60.8) color(gs12) clwidth(.8)) ///
					(function y=40, range(38 59.8) color(gs12) clwidth(.8)) ///
					(function y=39.5, range(38 59.5) color(gs12) clwidth(.8)) ///
					(function y=39, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38.5, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=38, range(39 59.3) color(gs12) clwidth(.8)) ///
					(function y=37.5, range(39 58.8) color(gs12) clwidth(.8)) ///
					(function y=37, range(40 58.3) color(gs12) clwidth(.8)) ///
					(function y=36.5, range(40 57.8) color(gs12) clwidth(.8)) ///
					(function y=36, range(41 56.8) color(gs12) clwidth(.8)) ///
					(function y=35.5, range(42 55) color(gs12) clwidth(.8)) ///
					(function y=35, range(42 55) color(gs12) clwidth(.8)) ///
					(scatter prop_y prop_x if treatment==3 & period==2 & proposer_dec==1 & prop_ucset==1, msymbol(D) mfcolor(gs8) msize(medsmall) mlcolor(black)) ///
					(scatter prop_y prop_x if treatment==3 & period==2 & proposer_dec==1 & prop_ucset==0, msymbol(S) mfcolor(gs14) msize(medsmall) mlcolor(black)) ///
					(scatteri 35 22 (6) "A", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 17 (12) "B", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 51 49 (12) "C", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 45 82 (12) "D", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					(scatteri 36 80 (6) "E", mcolor(black) msize(large) msymbol(o) mlabsize(medium) mlabcolor(black)) ///
					if treatment==3, ///     
						legend(size(vsmall) symxsize(4) order(1 "Uncovered Set"  34 "Proposal inside Uncovered Set"  35 "Proposal outside Uncovered Set") cols(3) rows(1) ring(0) bplacement(neast)) ///
						title(`"{bf:Treatment 3, Period 2}"', size(medium) color(black)) ///
						xtitle("") ///
						ytitle("") ///
						yscale(range(0 100)) ylabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						xscale(range(0 100)) xlabel(0 (10) 100, labsize(small) grid glwidth(thin)) ///    
						graphregion(color(white) margin(vsmall)) ///
						plotregion(lcolor(black) lwidth(medthick) margin(vsmall)) ///
						note("1 of 44 (2.3%) proposals inside" "Uncovered Set" "Average proposal: (42.2|46.1)", size(small)) ///
						xsize(4) ysize(4) ///
						name(Treatment3P2, replace)
				
			* COMBINE ALL GRAPHS	
				grc1leg Treatment1P1 Treatment1P2 Treatment2P1 Treatment2P2 Treatment3P1 Treatment3P2, ///
				graphregion(color(white) margin(vsmall)) ///
				graphregion(color(white) margin(vsmall)) ///
				cols(2) ///
				name(Figure_5, replace)
				graph display, ysize(14) xsize(7.8) scale(1) 

				graph close Treatment1P1 Treatment1P2 Treatment2P1 Treatment2P2 Treatment3P1 Treatment3P2 
		
		
		
********************************************************************************
** Supplementary Materials D - Outcomes of Group Decisions 						                
********************************************************************************				

	quietly do "Data_Reading_subjects.do"				
	
	sort treatment period groupid				
	list groupid final_x final_y final_ucset fin_points_a fin_points_b fin_points_c fin_points_d fin_points_e vote_a_adj vote_b_adj vote_c_adj vote_d_adj vote_e_adj time_in_stage2 num_ag if type==1
	

	
log close
